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Abstract 

We consider an elliptic self-adjoint first order pseudodifferential operator 
acting on columns of complex-valued half-densities over a connected com- 
pact manifold without boundary. The eigenvalues of the principal symbol 
are assumed to be simple but no assumptions are made on their sign, so 
the operator is not necessarily semi-bounded. We study the following ob- 
jects: the propagator (time-dependent operator which solves the Cauchy 
problem for the dynamic equation), the spectral function (sum of squares of 
Euclidean norms of eigenfunctions evaluated at a given point of the man- 
ifold, with summation carried out over all eigenvalues between zero and a 
positive A) and the counting function (number of eigenvalues between zero 
and a positive A). We derive explicit two-term asymptotic formulae for all 
three. For the propagator "asymptotic" is understood as asymptotic in terms 
of smoothness, whereas for the spectral and counting functions "asymptotic" 
is understood as asymptotic with respect to A — > +oo. 
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1 Main results 

The aim of the paper is to extend the classical results of [5] to systems. We are 
motivated by the observation that, to our knowledge, all previous publications on 
systems give formulae for the second asymptotic coefficient that are either incor- 
rect or incomplete (i.e. an algorithm for the calculation of the second asymptotic 
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coefficient rather than an actual formula). The appropriate bibliographic review 
is presented in Section [TT1 

Consider a first order classical pseudodifferential operator A acting on columns 
v = (vi . . . v m j of complex- valued half-densities over a connected compact n- 
dimensional manifold M without boundary. Throughout this paper we assume 
that m, n > 2 . 

We assume the coefficients of the operator A to be infinitely smooth. We also 
assume that the operator A is formally self-adjoint (symmetric): J M w*Avdx = 
J M (Aw)*v dx for all infinitely smooth v,w : M — > C m . Here and further on the 
superscript * in matrices, rows and columns indicates Hermitian conjugation in 
C m and dx := dx 1 . . . dx n , where x = (x 1 , . . . , x n ) are local coordinates on M. 

Let Ai(x,£) be the principal symbol of the operator A. Here £ = (£1,... , £ n ) 
is the variable dual to the position variable x; in physics literature the £ would 
be referred to as momentum. Our principal symbol A\ is an m x m Hermitian 
matrix-function on T'M := T*M \ {£ = 0}, i.e. on the cotangent bundle with the 
zero section removed. 

Let hfj' (x, £) be the eigenvalues of the principal symbol. We assume these 
eigenvalues to be nonzero (this is a version of the ellipticity condition) but do not 
make any assumptions on their sign. We also assume that the eigenvalues (x, £) 
are simple for all (x, £) £ T'M. The techniques developed in our paper do not work 
in the case when eigenvalues of the principal symbol have variable multiplicity, 
though they could probably be adapted to the case of constant multiplicity different 
from multiplicity 1. The use of the letter "h" for an eigenvalue of the principal 
symbol is motivated by the fact that later it will take on the role of a Hamiltonian, 
see formula (jl.f 31) . 

We enumerate the eigenvalues of the principal symbol hV'(x,£) in increasing 
order, using a positive index j = 1, . . . ,m + for positive h^'(x,^) and a negative 
index j = —1, . . . , — mT for negative h^'(x, £). Here m + is the number of positive 
eigenvalues of the principal symbol and mT is the number of negative ones. Of 
course, m + + mT = m. 

Under the above assumptions A is a self-adjoint operator, in the full functional 
analytic sense, in the Hilbert space L 2 (M; C™ 1 ) (Hilbert space of square integrable 
complex-valued column "functions") with domain H 1 (M;C m ) (Sobolev space of 
complex- valued column "functions" which are square integrable together with their 
first partial derivatives) and the spectrum of A is discrete. These facts are easily 
established by constructing the parametrix (approximate inverse) of the operator 
A + il. 

Let Afe and Vk — (vki(x) . . . Vk m {xj) be the eigenvalues and eigenfunctions 
of the operator A. The eigenvalues are enumerated in increasing order with 
account of multiplicity, using a positive index fc — f,2,... for positive and 
a nonpositive index k = 0,-1,-2,... for nonpositive Afc. If the operator A is 
bounded from below (i.e. if mT = 0) then the index fc runs from some integer 
value to +oo; if the operator A is bounded from above (i.e. if m + = 0) then the 
index k runs from — oo to some integer value; and if the operator A is unbounded 
from above and from below (i.e. if m + ^ and mT ^ 0) then the index k runs 
from — oo to +oo. 
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We will be studying the following three objects. 



Object 1. Our first object of study is the propagator, which is the one- 
parameter family of operators defined as 



U(t):=e- itA = Y,e- itXk v k (x) I [v k (y)]* ( • ) dy , (1.1) 
k Jm 



t e R. The propagator provides a solution to the Cauchy problem 

w \t=o = « 

for the dynamic equation 



(1.2) 



D t w + Aw = 0, (1.3) 

where D t := —id/dt. Namely, it is easy to see that if the column of half-densities 
v = v(x) is infinitely smooth, then, setting w := U(t) v, we get a time-dependent 
column of half-densities w(t, x) which is also infinitely smooth and which satisfies 
the equation (II .3[) and the initial condition (II. 2p . The use of the letter "IT" for the 
propagator is motivated by the fact that for each t the operator U (t) is unitary. 

Object 2. Our second object of study is the spectral function, which is the 
real density defined as 

e(X,x,x):= IM^II 2 - (1-4) 

0<A fc <A 

where ||i;fc(a;)|| 2 := [ v k(x)]*Vk(x) is the square of the Euclidean norm of the eigen- 
function i>k evaluated at the point x S M and A is a positive parameter (spectral 
parameter). 

Object 3. Our third and final object of study is the counting function 



N(\):= 1 = / e{X,x,x)dx. (1.5) 

0<A fc <A 



In other words, iV(A) is the number of eigenvalues Afe between zero and A. 

It is natural to ask the question: why, in defining the spectral function (|1.4[) and 
the counting function (jl.5l) , did we choose to perform summation over all positive 
eigenvalues up to a given positive A rather than over all negative eigenvalues up to 
a given negative A? There is no particular reason. One case reduces to the other 
by the change of operator A H> —A. This issue will be revisited in Section [TOT 

Further on we assume that m + > 0, i.e. that the operator A is unbounded 
from above. 



Our objectives are as follows. 



3 



Objective 1. We aim to construct the propagator (jl.ip explicitly in terms of 
oscillatory integrals, modulo an integral operator with an infinitely smooth, in the 
variables t, x and y, integral kernel. 

Objectives 2 and 3. We aim to derive, under appropriate assumptions on 
Hamiltonian trajectories, two-term asymptotics for the spectral function (|1.4[) and 
the counting function (jl.5[) . i.e. formulae of the type 

e(A, x, x) = a{x) A™ + b(x) A"" 1 + o(A™- 1 ), (1.6) 

N(X) = a\ n + bX"- 1 + o(A"- 1 ) (1.7) 

as A — > +oo. Obviously, here we expect the real constants a, b and real densities 
a{x), b{x) to be related in accordance with 

a = a(x) dx, (1-8) 

J M 

b= I b(x)dx. (1.9) 

It is well known that the above three objectives are closely related: if one 
achieves Objective 1, then Objectives 2 and 3 follow via Fourier Tauberian theo- 
rems ei nn m nzi • 

We are now in a position to state our main results. 

Result 1. We construct the propagator as a sum of m oscillatory integrals 

U ( t) ™*?~^ u V) ( j t)i (L10 ) 

3 

where the phase function of each oscillatory integral U^(t) is associated with 
the corresponding Hamiltonian h^' (x , £) . The symbol of the oscillatory inte- 
gral U^\i) is a complex- valued m x m matrix-function u^(t;y,rj), where y = 
(y 1 , . . . ,y n ) is the position of the source of the wave (i.e. this is the same y that 
appears in formula (ll.ip ) and r\ — . . . , ry„) is the corresponding dual variable 
(covector at the point y). When \rj\ —> +oo, the symbol admits an asymptotic 
expansion 

u^(t;y, V ) = 4 j) (t; y, V ) + (t;y, ?,) + ... (1.11) 

into components positively homogeneous in rj, with the subscript indicating degree 
of homogeneity. 

The formula for the principal symbol of the oscillatory integral (t) is known 
|16( [TTj and reads as follows: 

u^(t;y,r,) = [v^(x^(t;y,r,U^(t;y,r,))] [v^(y, V )}* 

xexp^| q^(x^(r;y,r,),^(T;y,r ] ))dT\ , (1.12) 
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where v^'(z,Q is the normalised eigenvector of the principal symbol A\{z, £) 
corresponding to the eigenvalue (Hamiltonian) hP> (z, £), [x^> (t; y, Tj),^) (f- ^ 77)) 
is the Hamiltonian trajectory originating from the point (y, 77), i.e. solution of the 
system of ordinary differential equations (the dot denotes differentiation in t) 

A® = h { p(x^,^), = -^(i^,^ 1 ) (1-13) 

subject to the initial condition (x (j) , \ t=Q = (y,T)), <? (j) : T'M -> R is the 
function 

q <J) ■- [ v (i)]*A suh v^ - ^{[vW]*,Ai-hW,vW}-i[ v W]*{ v W,hV>} (1.14) 
and 

A suh (z, C) := A (z, C) + \{A x ) zHa (z, C) (1.15) 

is the subprincipal symbol of the operator A, with the subscripts z a and C, a in- 
dicating partial derivatives and the repeated index a indicating summation over 
a = l,...,n. Curly brackets in formula (|1.14|) denote the Poisson bracket on 
matrix-functions 

{PR} := P z «R Qa - P Ca R z « (1.16) 
and its further generalisation 

{P,Q,R}:=P z «QR Ca -P ia QR z «. (1.17) 

As the derivation of formula (|1.12j) was previously performed only in theses 
|16[ 111), we repeat it in Sections [2] and [3] of our paper. Our derivation differs 
slightly from that in [16] and |11) . 

Formula (jTTT3J) is invariant under changes of local coordinates on the manifold 
M, i.e. elements of the mxm matrix-function Uq (t; y, rj) are scalars on R x T'M. 
Moreover, formula (|1.12|) is invariant under the transformation of the eigenvector 
of the principal symbol 

V U) e ** W) u 0) ) (1.18) 

where 

W : T'M -> K (1.19) 

is an arbitrary smooth function. When some quantity is defined up to the action 
of a certain transformation, theoretical physicists refer to such a transformation as 
a gauge transformation. We follow this tradition. Note that our particular gauge 
transformation (|1.18p . (|1.19p is quite common in quantum mechanics: when <pU) 
is a function of the position variable x only (i.e. when <j)^' : M — > R) this gauge 
transformation is associated with electromagnetism. 

Both Y. Safarov [TC] and W.J. Nicoll [TT] assumed that the operator A is semi- 
bounded from below but this assumption is not essential and their formula (|1.12D 
remains true in the more general case that we are dealing with. 

However, knowing the principal symbol (|1.12[) of the oscillatory integral U^\t) 
is not enough if one wants to derive two-term asymptotics fj 1 . 6|) and (|1.7jl . One 
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needs information about u_{(t; y, rf), the component of the symbol of the oscilla- 
tory integral U^\t) which is positively homogeneous in 77 of degree -1, see formula 

(jl.lljl . but here the problem is that y, if) is not a true invariant in the sense 

that it depends on the choice of phase function in the oscillatory integral. We 
overcome this difficulty by observing that U^> (0) is a pseudodifferential operator, 
hence, it has a well-defined subprincipal symbol \U^' (0)] su b- We prove that 

tr[C/«(0)] sub = -i{[^)]*,«^} (1.20) 

and subsequently show that information contained in formulae (|1.12jl and (11.201) 
is sufficient for the derivation of two-term asymptotics (|1.6p and (|1.7[1 . 

Note that the RHS of formula (|1.20p is invariant under the gauge transforma- 
tion (fTTTHl) . (091) . 

Formula (|1.20|1 plays a central role in our paper. Sections [2] and [3] provide 
auxiliary material needed for the proof of formula (|1.20| , whereas the actual proof 
of formula (|1 .20|) is given in Section |4] 

Let us elaborate briefly on the geometric meaning of the RHS of (|1.20p (a more 
detailed exposition is presented in Section [5J . The eigenvector of the principal 
symbol is defined up to a gauge transformation (|1.18|l . (|1.19[) so it is natural to 
introduce a U(l) connection on T'M as follows: when parallel transporting an 
eigenvector of the principal symbol along a curve in T'M we require that the 
derivative of the eigenvector along the curve be orthogonal to the eigenvector 
itself. This is equivalent to the introduction of an (intrinsic) electromagnetic field 
on T'M , with the 2n-component real quantity 

t([v^]*v { j] , [v&]*vg) (1.21) 

playing the role of the electromagnetic covector potential. Our quantity (jl.21|) is 
a 1-form on T'M, rather than on M itself as is the case in "traditional" electro- 
magnetism. The above U(l) connection generates curvature which is a 2-form on 
T'M , an analogue of the electromagnetic tensor. Out of this curvature 2-form one 
can construct, by contraction of indices, a real scalar. This scalar curvature is the 
expression appearing in the RHS of formula (|1.20|) . 

Observe now that U^(0) is the identity operator on half-densities. The 
subprincipal symbol of the identity operator is zero, so formula (|1.20p implies 

J2{[v U) ]*,v^} = 0. (1.22) 

3 

One can check the identity (|1.22l) directly, without constructing the oscillatory 
integrals U^\t): it follows from the fact that the v^'{x,^) form an orthonormal 
basis, see end of Section [5] for details. We mentioned the identity (jl.22p in order 
to highlight, once again, the fact that the curvature effects we have identified are 
specific to systems and do not have an analogue in the scalar case. 
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Results 2 and 3. We prove, under appropriate assumptions on Hamiltonian 
trajectories (see Theorems 18.31 and 18. 4p . asymptotic formulae (|1.6[) and (|1.7p with 




(1.23) 



- i{[ v W]*,A 1 -h&,vV ) }+-^hW{[vW] m ,vW}yx i €)aZ, (1.24) 

and a and b expressed via the above densities (|1.23l) and (|1.24[) as (|1.8[) and (|1.9p . 
In (|I~23)) and (fl~24| ^ is shorthand for <Z£ := (2tt)-" ^ = (27r) -n d& . . . d£ n , 
and the Poisson bracket on matrix-functions { • , • } and its further generalisation 
{ ■ , • , • } are defined by formulae (| 1 . 16[) and (|1.17p respectively. 

To our knowledge, formula (| 1 . 24|) is a new result. Note that in [16] this for- 
mula (more precisely, its integrated over M version (|1.9p ) was written incorrectly, 
without the curvature terms — -^4- J h^{[v^]*, v^'}. See also Section HT1 where 
we give a more detailed bibliographic review. 

It is easy to see that the right-hand sides of (| 1 . 23|) and f| 1 . 24|) behave as densities 
under changes of local coordinates on the manifold M and that these expressions 
are invariant under gauge transformations (|1.18p , (|1.19p of the eigenvectors of the 
principal symbol. Moreover, the right-hand sides of (I1.23P and (jl.24p are unitarily 
invariant, i.e. invariant under the transformation of the operator 

A^RAR*, (1.25) 

where 

R: M -> U(m) (1.26) 

is an arbitrary smooth unitary matrix-function. The fact that the RHS of (|1.24p is 
unitarily invariant is non-trivial: the appropriate calculations are presented in Sec- 
tion[9l The observation that without the curvature terms — ^zj f h^'{[v^']*,v^'} 
(as in |16| ) the RHS of (jl.24p is not unitarily invariant was a major motivating 
factor in the writing of this paper. 

Formula (fl~24|) is the main result of our paper. Note that even though the 
two-term asymptotic expansion (|1.6p holds only under certain assumptions on 
Hamiltonian trajectories (loops), the second asymptotic coefficient f| 1 . 24f) is, in 
itself, well-defined irrespective of how many loops we have. If one wishes to refor- 
mulate the asymptotic expansion (jl.6p in such a way that it remains valid without 
assumptions on the number of loops, this can easily be achieved, say, by taking a 
convolution with a function from Schwartz space <S(R): see Theorem 17. II 
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2 Algorithm for the construction of the propagator 

We construct the propagator as a sum of m oscillatory integrals (| 1 . 10[) where each 
integral is of the form 

uV)(t) = J e^^^u^H^y^^i^x^^d^^x^^i^dydr,. (2.1) 

Here we use notation from the book [T3J, only adapted to systems. Namely, the 
expressions appearing in formula (|2.1[) have the following meaning. 

• The function ipV> is a phase function, i.e. a function K x M x T'M — > C 
positively homogeneous in rj of degree 1 and satisfying the conditions 

<p® (t, x; y, v) = (x~ x^ (t; y, rijf (t; y, r,) +0(\x- x& (t; y, r,) | 2 ), (2.2) 

Im tp {j) (t, x; y, rj) > 0, (2.3) 

det cp% fj (t,xW(t;y ) riy,y,r ] )^0. (2.4) 

Recall that according to Corollary 2.4.5 from [IS] we are guaranteed to have 
(|2.4p if we choose a phase function 

<p® (t, x; y, V) = (x~ (t; y, r,)) a & ft y, v) 

+ \c% (t; V, V) (x - (t; y, V )T (x - (t; y, r,)f 

+ 0(\x-x^(t-y, V )\ 3 ) (2.5) 

with complex-valued symmetric matrix-function C^g satisfying the strict 
inequality Im > (our original requirement (|2.3[) implies only the non- 
strict inequality ImC® > 0). Note that even though the matrix- function 
C^p is not a tensor, the inequalities Im C^' > and Im > are invariant 
under transformations of local coordinates x; see Remark 2.4.9 in [18] for 
details. 

• The quantity u^> is the symbol of our oscillatory integral, i.e. a complex- 
valued m x m matrix-function K x T'M — > C m " which admits the asymptotic 
expansion p. lip . The symbol is the unknown quantity in our construction. 

• The quantity d^j) is defined in accordance with formula (2.2.4) from [THJ as 

d^&x-^n) := (detV&jV* = \detS m \ 1/2 e i - s(det2 ^ )/i . (2.6) 

Note that in view of (|2.4p our eLyj is well-defined and smooth for x close 
to x^i (t;y,r)). It is known |18| that under coordinate transformations gLo) 
behaves as a half-density in x and as a half-density to the power —1 in y. 

In formula (|2.6p we wrote (det 2 ^^.) 1 / 4 rather than (det V 3 !""^) 1 ^ 2 i n order 
to make this expression truly invariant under coordinate transformations. 
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Recall that local coordinates x and y are chosen independently and that rj 
is a covector based at the point y. Consequently, det changes sign 

under inversion of one of the local coordinates x a , a — l,.,.,n, or y@, 
(3 = 1, . . . , n, whereas det 2 ^}^ retains sign under inversion. 

The choice of (smooth) branch of arg^et 2 ^!^ ) is assumed to be fixed. 
Thus, for a given phase function ip^> formula (|2.6p defines the quantity 
d v (j) uniquely up to a factor e lk7T / 2 : k = 0,1,2,3. Observe now that if 
we set t = and choose the same local coordinates for x and y, we get 
V 3 £""77 /3 (0i 2/^ 7 ?) = I- This implies that we can fully specify the choice of 
branch of arg(det 2 y^L JJ ) by requiring that d^u) (0, y,y,T}) = 1. 
The purpose of the introduction of the factor d v u) m (|2.ip is twofold. 

(a) It ensures that the symbol is a function on E x T'M in the full 
differential geometric sense of the word, i.e. that it is invariant under 
transformations of local coordinates x and y. 

(b) It ensures that the principal symbol Uq does not depend on the choice 
of phase function ip^\ See Remark 2.2.8 in [TS] for more details. 

• The quantity i s a smooth cut-off function RxA/x T'M -» M satisfying 
the following conditions. 

(a) ^ {t, x; y, V ) =0 on the set {(t,x;y,r,): \h^ (y, rf)\ < 1/2}. 

(b) ?W (t, x; y, rj) = 1 on the intersection of a small conic neighbourhood of 
the set 

{(t,x;y, V ): x = x^(t;y,r})} (2.7) 
with the set {(t, x;y, 77) : \h^(y,rj)\ > 1}. 

(c) ^(t,x;y,\r)) =^(t,x;y,rf) for |/i (j %,r7)| > 1, A > 1. 

• It is known (see Section 2.3 in [TB] for details) that Hamiltonian trajectories 
generated by a Hamiltonian (x, £) positively homogeneous in £ of degree 1 
satisfy the identity 

{x U) r ^U) =0i (2 . 8) 
where (4 j) ) Q/3 := d(x^) a /dr]p. Formulae JMH and JMJ) imply 

cp^(t,x^(t;y,r,);y,r 1 )=0. (2.9) 

This allows us to apply the stationary phase method in the neighbourhood 
of the set (|2.7[) and disregard what happens away from it. 



Our task now is to construct the symbols Uq (t; y, rf), j — 1, . . . , to, so that our 
oscillatory integrals U^'(t), j = 1, . . . , to, satisfy the dynamic equations 

(D t +A(x,D x ))U^(t) mod = C °° (2.10) 
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and initial condition 

J2uW(0) mod =°° I, (2.11) 

3 

where / is the identity operator on half-densities; compare with formulae (jl.3p . 
(|1.2p and Ql.lOp . Note that the pseudodifferential operator A in formula (|2.10p acts 
on the oscillatory integral U(t) in the variable x; say, if A is a differential operator 
this means that in order to evaluate AU^\t) one has to perform the appropriate 
differentiations of the oscillatory integral (|2.ip in the variable x. Following the 
conventions of Section 3.3 of [T8], we emphasise the fact that the pseudodifferential 
operator A in formula (|2.10p acts on the oscillatory integral U(t) in the variable x 
by writing this pseudodifferential operator as A(x,D x ), where D x a := —id/dx a . 
We examine first the dynamic equation (|2.10p . We have 

(D t + A(x,D x ))U^(t)=F^(t), 

where F^> (t) is the oscillatory integral 

( t ) = J e^ (j) / W (t, x- y, ri) (i, x; y, r?) d y w (t, re; y, rj) ( • ) dy d V 

whose matrix-valued amplitude f^> is given by the formula 

f U) = DtU U) + (^0) + (d (pW )- 1 (D t d (pW ) + u«, (2.12) 

where the matrix-function sw (t, x; y, 77) is defined as 

S U) = e -^ U) (d vU) r 1 A(x,D x )(e^ U) d v(j) ). (2.13) 

Theorem 18.1 from |19| gives us the following explicit asymptotic (in inverse 
powers of 77) formula for the matrix-function (|2.13p : 

= (d^r'Y.^^^^i.Dfx^)]^ , (2.14) 

where 

X W (t,z,x;y,T,) = e i ^ t ^^d v u)(t,z ] y,T,), (2.15) 

Tp (j) (t,z,x;y,v) = <P ij) (t,z;y,r)) - ip U) (t, x; y, rj) - ^ (i, x; y, 77) (z - x) 13 . (2.16) 
In formula (|2.14jl 

• a := (ai,...,a n ) is a multi-index (note the bold font which we use to 
distinguish multi-indices and individual indices), ot\ := a\\ ■ ■ ■ a n \ , D" := 

D"l ■ ■ ■ D™£, D z p := —id/dzP, 

• A(x, £) is the full symbol of the pseudodifferential operator A written in local 
coordinates x, 

. A^(x,0 := d?A(x,$, 3f := S£ ■ ■ ■ and d u := d/d£p . 
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When \rj\ — > +00 the matrix- valued amplitude f^\t, x; y, 77) defined by formula 
(12. 12|) admits an asymptotic expansion 

(t, x; y, n) = f[ j) (t, x- y, r,) + (t, x; y, 77) + /«> (t, x;y,rf) + ... (2.17) 

into components positively homogeneous in 77, with the subscript indicating degree 
of homogeneity. Note the following differences between formulae (jl.lip and (|2.17p . 

• The leading term in (|2.17[) has degree of homogeneity 1, rather than as in 
(II. lip . In fact, the leading term in (I2.17P can be easily written out explicitly 

f[ j) {t, x; y, 77) = fa® (t, x- y, rj) + A x (x, <pV) ( t) X] y , u « ( t; ^ v ) , (2.18) 

where Ai(x,^) is the (matrix- valued) principal symbol of the pseudodiffer- 
ential operator A. 

• Unlike the symbol u^' (t; y , 77) , the amplitude (t,x;y,rj) depends on x. 

We now need to exclude the dependence on x from the amplitude (t, x; y, rf). 
This can be done by means of the algorithm described in subsection 2.7.3 of [18J. 
We outline this algorithm below. 

Working in local coordinates, define the matrix-function ipxr/ in accordance 

with ((pxn)a^ '■= fx"riB an< i then define its inverse (y4"§) _1 from the identity 

{ L P^)a^[{ l Pxri)^ 1 ]i} 1 '■— 5a 1 . Define the "scalar" first order linear differential oper- 
ators 

4 J) V(0/aA a=l,...,n. (2.19) 

Note that the coefficients of these differential operators are functions of the position 
variable x and the dual variable £. It is known, see part 2 of Appendix E in |18) . 
that the operators (|2.19p commute: La — L^La, ot,P — 1, . . . ,n. 

Denote L { £ := (L«)^ • • • (L$) a «, {-^T ■= • • • (-^) a ", 

and, given an r G N, define the "scalar" linear differential operator 

where \a\ := ol\ + . . . + a n and the repeated index /3 indicates summation over 
/3 = l,...,n. 

Recall Definition 2.7.8 from [18]: the linear operator L is said to be positively 
homogeneous in 77 of degree p £ R if for any q G K and any function / positively 
homogeneous in 77 of degree q the function Lf is positively homogeneous in 77 of 
degree p + q. It is easy to see that the operator (I2.20p is positively homogeneous 

(i) 

in 77 of degree —1 and the first subscript in r emphasises this fact. 
Let &^ be the (linear) operator of restriction to x — (t;y, 77), 

'So''' : = i')\x = xU)(t-y,T 1 ) ! ( 2 - 21 ) 
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and let 



e U) r ■.= e ( j) (^Lr (2.22) 



( i) 

for r = 1,2,.... Observe that our linear operators &l r , r — 0, 1, 2, . . ., are pos- 
itively homogeneous in r\ of degree — r. This observation allows us to define the 
linear operator 

+ OG 

6« := &U l > ( 2 - 2 3) 

r=0 

where the series is understood as an asymptotic series in inverse powers of r/. 

According to subsection 2.7.3 of [IS], the dynamic equation (|2.10p can now be 
rewritten in the equivalent form 

6W)/0')= o, (2.24) 

where the equality is understood in the asymptotic sense, as an asymptotic expan- 
sion in inverse powers of rj. Recall that the matrix- valued amplitude /00 (t, x\ y, n) 
appearing in (j!T2"4"| is defined by formulae ([2~T2]) - ([2TT51) . 

Substituting (|2.23[) and f|2 . 1 7|) into (|2.24p we obtain a hierarchy of equations 

e^f^ = o, (2.25) 
6 W/(i) + 6 U)/W =0) (226) 

6^+6^+60^ = 0, 



positively homogeneous in 7/ of degree 1, 0, —1, .... These are the transport equa- 
tions for the determination of the unknown homogeneous components (t; y, 77), 
u_i(t; y, 77), u^ 2 {t\ Hi 7 ]): ■ ■ •! of the symbol of the oscillatory integral (|2.1D . 

Let us now examine the initial condition (|2.11[) . Each operator U^'(Q) is a 
pseudodifferential operator, only written in a slightly nonstandard form. The 
issues here are as follows. 

• We use the invariantly defined phase function (^00 (0, x\ y, rf) — (x — y) a r\ a + 
OQx — y\ 2 ) rather than the linear phase function [x — y) a r\ a written in local 
coordinates. 

• When defining the (full) symbol of the operator U^\t) we excluded the 
variable x from the amplitude rather than the variable y. Note that when 
dealing with pseudodifferential operators it is customary to exclude the vari- 
able y from the amplitude; exclusion of the variable x gives the dual symbol 
of a pseudodifferential operator, see subsection 2.1.3 in [T8]. Thus, at t = 0, 
our symbol it 00 (0; y, 77) resembles the dual symbol of a pseudodifferential 
operator rather than the "normal" symbol. 



We have the extra factor d v u) (0, x; y, rf) in our representation of the operator 
t/00(0) as an oscillatory integral. 
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The (full) dual symbol of the pseudodifferential operator LTW (0) can be calcu- 
lated in local coordinates in accordance with the following formula which addresses 
the issues highlighted above: 

E { ~^- ( D " d " uU) (o; * v) eiuW {x ' ,v ' n) d ^ (°> x ' y> ♦/)) L=, . ( 2 - 27 ) 

where a;^ (x; y, 77) = <^( J ) (0, x; y, 77) — (x — y)' 9 77^ . Formula (|2.27p is a version of 
the formula from subsection 2.1.3 of [TB], only with the extra factor (— l)\ a K The 
latter is needed because we are writing down the dual symbol of the pseudodiffer- 
ential operator U(°'(Q) (no dependence on x) rather than its "normal" symbol (no 
dependence on y). 

The initial condition p. lip can now be rewritten in explicit form as 

E E (jD " d " uU) (0; y ' v) e " <3> {x ' y ' v) dvU) (0 ' x; Vi v) ^=y = 1 ' (2 - 28) 

j a 

where / is the m x m identity matrix. Condition (|2.28p can be decomposed into 
components positively homogeneous in 77 of degree 0, — 1, —2, . . ., giving us a hierar- 
chy of initial conditions. The leading (of degree of homogeneity 0) initial condition 
reads 

J2^\0;y,r))=I, (2.29) 

i 

whereas lower order initial conditions are more complicated and depend on the 
choice of our phase functions ipV'. 

3 Leading transport equations 

Formulae (j2~2"Tj) . (|2~T5j) . (l2T2"l) . (TElS)) and the identity £ a h^(x,£) = h^(x,£) 
(consequence of the fact that hV' (x, £) is positively homogeneous in £ of degree 1) 
give us the following explicit representation for the leading transport equation 

[A, (z« (t; y, 77), (t; y, 77)) -h& (t; y, 77), ^ (t; 2/, »?))] 4'' 3 (*5 V, l) = «• (3.1) 

Here, of course, h^' (x^'(t;y,rf),^'(t;y,rf)) = h^\y,rf). 
Equation Q3.ip implies that 

v$ (t; y, 77) = u« (ajW (i; y, 77), (t; y, 77)) (t; y, V )f , (3.2) 

where (z, is the normalised eigenvector of the principal symbol A\{z,C) 
corresponding to the eigenvalue h^'(z, £) and w<- j '> : M x T'M -> C m is a column- 
function, positively homogeneous in 77 of degree 0, that remains to be found. 
Formulae (|2.29p and (|3.2p imply the following initial condition for the unknown 

column- function w^' : 

v>W(0;y,n) = v&(y,Ti). (3.3) 
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We now consider the next transport equation in our hierarchy, equation (|2.26p . 
We will write down the two terms appearing in (|2 . 26[) separately. 

In view of formulae (|2~Tgj) and (|2^D| - (|2~2^1) . the first term in (l2~2l)l) reads 



(3.4) 



where we dropped, for the sake of brevity, the arguments (<; y, rj) in Uq and 
x^\ and the arguments {t,x;y,rj) in (p^\ <fx\ frP an d eLo) . Recall that the 
differential operators La are defined in accordance with formula (|2.19l) and the 
coefficients of these operators depend on (t, x; y, rf). 

In view of formulae (|2TT2|) - (|2TT7|1 and (12T2T1) . the second term in (l2~2l))) reads 



(=0) f U) 



(i) 



(d„ )) _1 (A + (-Ai)^^-)^ + A 



,00 



x=xU) 

+ [A 1 -h®]uV\, (3.5) 



where 



(3.6) 



is the matrix-function from (|2.5[) . In formulae (13.51) and (|3.6p we dropped, for the 
sake of brevity, the arguments (t;y,rj) in it^, w'fj, C^2 and x^\ the arguments 
(x^\t;y, r]),^(t;y, ij)) in A , A x , (Ai)^ a , (Ai)^ p and /iW, and the arguments 

(t, x; y, T]) in d^j) and (/^j^ . 

Looking at (|3.4p and (|3 . 5|) we see that the transport equation (|2.26p has a 
complicated structure. Hence, in this section we choose not to perform the analysis 
of the full equation ()2.26|) and analyse only one particular subequation of this 
equation. Namely, observe that equation (|2.26ll is equivalent to m subequations 



[6 uyu) 



6 W /o W ] 



0, 



1*3, 



(3.7) 
(3.8) 



where we dropped, for the sake of brevity, the arguments (x^ (t; y, 77), ft?') (t; y, 77)) 
in lyv)] and [v^] . In the remainder of this section we analyse (sub)equation 
(|3~7f only. 

Equation ()3.7p is simpler than each of the m — 1 equations (|3 . 8[) for the following 
two reasons. 

• Firstly, the term [Ai — h^]u^\ from (13. 5p vanishes after multiplication by 
[oW] * from the left. Hence, equation (|3.7p does not contain u"l\. 
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Secondly, if we substitute p. 21) into (|3.7[) . then the term with 
vanishes. This follows from the fact that the scalar function 

has a second order zero, in the variable x, at x = x^(t;y,rf). Indeed, we 
have 
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9 ^b (5 T(*'r+^(^i j) ))« w 

_ (#+A 1 (x,^))) ;e 



= - Cgftg + + cg(A l)? > W 

= [«W]'(Ax)..i,W) -fcW +C<J([«W]*(A 1 ) e ,t,W -*g>) =0, 

where in the last two lines we dropped, for the sake of brevity, the argu- 
ments (x^(t',y,ri),^\t',y,r))) in (A^c, (Ai)^, h^l, h^, and the argu- 
ment (t; y, rf) in cfl (the latter is the matrix-function from formulae (|2.5|1 
and (|3.6p ). Throughout the above argument we used the fact that our \v^] * 
and do not depend on x: their argument is (x^'(t;y, rf),^' (t;y, rf)). 



Substituting ([23), (GE3) and (g^J) into we get 

(A+P (j) (t;y,'7))M j) (i;y^)] T = o. 

where 



(3.9) 



p (3) = i [ v 00] 



^o-icAO^cg)^). (3.10) 



Note that the ordinary differential operator in the LHS of formula (|3.9p is a 
scalar one, i.e. it does not mix up the different components of the column-function 
w^\t] y, rf). The solution of the ordinary differential equation (|3.9p subject to the 
initial condition (|3.3p is 

(t;y f r)) = vW(y,rj)exi>(-ij^ p (l} (r;y,r?) dr^ . (3.11) 



Comparing formulae (|3.2p , (|3.1ip with formula (|1.12p we see that in order to prove 
the latter we need only to establish the scalar identity 



p {J Ht;y,r ] ) = q^(x^(t;y,r 1 ),^(t;y,r 1 )) : 



(3.12) 
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where qv> is the function (|1.14[) . In view of the definitions of the quantities 
and q(i\ see formulae p.lOp and (|1.14[) , and the definition of the subprincipal 
symbol (|1.15p . proving the identity (|3.12p reduces to proving the identity 



-2 [v^(x^,C U) )}* 



d 
drjp 



+ [^(a-CA ,£(*))]* ((Ax),^ + (A^CS) v<*HxW,tW). (3.13) 



Note that the expressions in the LHS and RHS of (|3.13p have different structure. 
The LHS of ([3~T3| is the generalised Poisson bracket {[v^]*,Ai - h^\v^}, see 
(fTTT7| . evaluated at z = x^{t;y,rf), ( = ^(t;y,rj), whereas the RHS of (|5?T5) | 
involves partial derivatives (in rf) of [x^' (t; y : rf),£^' (f- y ; 77)) (Chain Rule). In 
writing (13 . 1 3|) we also dropped, for the sake of brevity, the arguments (t, x; y, rf) 
in ipt , <fx j iptj , d v u) and the coefficients of the differential operators La an d 
L^j\ the arguments (x^\ in /il^, (^4i)a«^ a and (Ai)^^, and the arguments 
(t;y,r,) in and eg 



Before performing the calculations that will establish the identity (|3 . 1 3|) we 
make several observations that will allow us to simplify these calculations consid- 
erably. 

Firstly, our function pW) (t- y y rf) does not depend on the choice of the phase 
function <pV' (t,x;y,rf). Indeed, if p^(t; y, rf) did depend on the choice of phase 
function, then, in view of formulae p.2p and p. lip the principal symbol of our 
oscillatory integral U^\t) would depend on the choice of phase function, which 
would contradict Theorem 2.7.11 from [T8]. Here we use the fact that operators 
(t) with different j cannot compensate each other to give an integral operator 
whose integral kernel is infinitely smooth in t, x and y because all our U^\t) 
oscillate in t in a different way: tp^' (t,x^(t;y,r));y,rf) = —h^(y,rf) and we 
assumed the eigenvalues h^\y,rf) of our principal symbol Ax(jj, rf) to be simple. 

Secondly, the arguments (free variables) in (|3.13p are (t]y,rj). We fix an ar- 
bitrary point (t; y, fj) G K x T'M and prove formula (|3.13|) at this point. Put 
(C^)a /3 : = d(£^)a/dvp- According to Lemma 2.3.2 from [18] there exists a local 
coordinate system x such that det(£^) Q ^ ^ 0. This opens the way to the use of 
the linear phase function 



<p& (t, x- y, rf) = {x- x& (t; y, r,)) a & (t; y, r,) 



(3.14) 



which will simplify calculations to a great extent. Moreover, we can choose a local 
coordinate system y such that 



(ti j) )J(t-,y,rj) = Sj 



(3.15) 



which will simplify calculations even further. 
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The calculations we are about to perform will make use of the symmetry 

o#>f°(#v = o^r^v ( 3 - 16 ) 

which is an immediate consequence of formula (|2.8[) . Formula (|3.16[) appears as 
formula (2.3.3) in [18] and the accompanying text explains its geometric meaning. 
Note that at the point (t; y, fj) formula (|3.16ll takes the especially simple form 

{x^{t ] y,f,) = {xfY a {t-,y,^ (3.17) 

Our calculations will also involve the quantity <p^ave Vi v) where x := 
xW(t;y,fj). Formulae (j3~Hj) . (l2~g)) . (13~T5)) and ([337]) imply 

(3.18) 



^(t,x;y,ry) = -(xWr(t;y,ry) 



U)\<*Pt 



Further on we denote | := £w) (£; ^). 

With account of all the simplifications listed above, we can rewrite formula 
(|3.13l) . which is the identity that we are proving, as 



{[v^]*,A!-h U) ,v^}(x,i) 

d 2 



-{A 1 (x,^)-h^(y, V ) 



(x\?>y 



■ dx a drj a 

— [x — x 

^^(x^-h^v) 

— [x — x 



where ««) = = x^ (t; y , fj) , (AiJ.-c - (^i)*<*£ Q (s, I), fc^. 



(s,7j)=(2,7j) 

(3.19) 

0) _ 



^ = ^(5,0. = * W) and = €«)(t;y,»j). 



Note that the last two terms in the RHS of p,19|) originate from the term with 
d v (j) in p,13[) : we used the fact that d^u) does not depend on x and that 

K^r^t^lL ^ a-~~^ = -~(h% +h U l Jxi j) ) a P). (3.20) 

Formula (|3 . 20[) is a special case of formula (3.3.21) from |18j . 

Note also that the term —h^(y, if) appearing (twice) in the RHS of (13. 19)) will 



vanish after being acted upon with the differential operators qJLq^ 
because it does not depend on x. 
We have 



and 



d 2 

dx a dx@ 



teCi)i 



d 2 



dx a drj a 



{x,rf) = {x,fi) 



[^^(ii^^-^-^C^)^ 

+ [«»]*((i 1 ) x «-^2)(^+f;«(^)^), (3.21) 
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teCi)] 



d 2 
dx a dxP 



' (x,ri) = (x,fj) 

= [v^nA 1 ) xax ^\ (3.22) 



where (Ix)^ = (^(z^), j$ = /£2(z,f), gg> = vgJfoO and v$ = 
vjp(x,i). We also have 

[s w ]*((4)x« - ^2)g« + [g«r((^ - 

= ^-[g«]*(ixWgW. (3.23) 
Using formulae (13.231) and (|3.17p we can rewrite formula (|3.2ip as 



^-(^(x,^) - (x - xWyh$(xU),£W)) v^(x^,C U) ) 



U[v^r(A 1 ) xax ,v^ + ~h { X li ){xf)^. (3.24) 



Substituting flgJH) and (pT2"2"j) into ([3"T§|) we see that all the terms with {x^ ] ) al3 
cancel out and we get 



ft«)g«. 



(3.25) 



Thus, the proof of the identity p,13[) has been reduced to the proof of the iden- 
tity (pT2"5j) . 

Observe now that formula (|3.25|l no longer has Hamiltonian trajectories present 
in it. This means that we can drop all the tildes and rewrite Q3.25P as 



{[v&]*,Ai -h (j \v {l) } 



[v^]*(A 1 -h U) ) 



p) _ 2[v {j) ]*(Ai - h {3) ) a v[ j) , (3.26) 



where the arguments are [x, £). We no longer need to restrict our consideration to 
the particular point (x, £) = (x, £): if we prove (|3.26p for an arbitrary (x, £) G T'M 
we will prove it for a particular (i, £) £ T'M. 

The proof of the identity (|3.26p is straightforward. We note that 

[v {]) ]*{A l ~h^) x . ia v { ^ = 

- [v if) ]*(A 1 - h^) x «v£ - [v U} ]*(A 1 - h^) (a v^ (3.27) 
and substituting (I3.27P into (I3.26P reduce the latter to the form 

{[v^]\A 1 -h {j \v^} = 
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But 

(A! - hW) aa = -[v^ViA, ~ ftW)), (3.29) 

[« w ]*(Aa - ft^) fa = ~ hU) )- ( 3 - 30 ) 

Substituting (ET2"9"j) and $£3$ into (EQgf we get 

{[v&]*,Ai - hW )V V)} = [vi J }]*(A 1 - h&)vg - [v^]*(Ai - hV>)v$> 
which agrees with the definition of the generalised Poisson bracket (JTTT7J) . 



4 Proof of formula ( 11.201 ) 



In this section we prove formula (| 1.200 . Our approach is as follows. 
We write down explicitly the transport equations Q3.8[l at t = 0, i.e. 

[v d)y [6 0] / p-) + 6 a) / 0-)]|^ = 0j w . (4il) 

We use the same local coordinates for x and y and we assume all our phase 
functions to be linear, i.e. we assume that for each j we have (|3. 14[) . Using linear 
phase functions is justified for small t because we have )a^(0; Mi v) — and, 
hence, det V'i"^ x l V: v) f° r small t. Writing down equations (|4.1[) for linear 
phase functions is much easier than for general phase functions (|2.2p . 

Using linear phase functions has the additional advantage that the initial con- 
dition (|2.28p simplifies and reads now Y^j u V'(0; y, if) = I. In view of (jl.lip . this 
implies, in particular, that 

5>^(0)=0. (4-2) 

3 

Here and further on in this section we drop, for the sake of brevity, the arguments 

(y, v) in u -\- 

Of course, the formula we are proving, formula (|1.20p . does not depend on our 
choice of phase functions. It is just easier to carry out calculations for linear phase 
functions. 

We will show that (|4.ip is a system of complex linear algebraic equations for 
the unknowns u_\(0). The total number of equations (|4.ip is m 2 — m. However, 
for each j and I the LHS of (|4.1I) is a row of m elements, so (|4.ip is, effectively, a 
system of m(m 2 — m) scalar equations. 

Equation (|4.2p is a single matrix equation, so it is, effectively, a system of m 2 
scalar equations. 

Consequently, the system (|4.ip . (|4.2p is, effectively, a system of m 3 scalar equa- 
tions. This is exactly the number of unknown scalar elements in the m matrices 
u^(0). 

In the remainder of this section we write down explicitly the LHS of (|4.ip and 
solve the linear algebraic system (|4.ip . (14. 2p for the unknowns u^\(0). This will 
allow us to prove formula (|1.20p . 
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Before starting explicit calculations we observe that equations (|4.ip can be 
equivalently rewritten as 

P (i) [ 6 «)/W + 6 «)/«)]| t=o = 0, Ijtj, (4.3) 

where PW :— [yW (y, rj)} [v^ (y, rj)}* is the orthogonal projection onto the eigenspace 
corresponding to the (normalised) eigenvector «W (y, rj) of the principal symbol. 
We will deal with (|4.3p rather than with (14. ip . This is simply a matter of conve- 



4.1 Part 1 of the proof of formula ( I1.20B 

Our task in this subsection is to calculate the LHS of (|4.3[) . In our calculations 
we use the explicit formula (jl.121) for the principal symbol Ug\t;y,r)) which was 
proved in Section [21 

At t = formula (l3~4l reads 



<9 2 



(Ax (x, „) - h& (y,ri)-(x- V yh$ (y, V ))P U) (y, v) 



x=y 



which gives us 



U) 



(4.4) 



In the latter formula we dropped, for the sake of brevity, the arguments (y,r]). 
At t = formula (13. 5p reads 



CO 



o 



t=o 



-i{vV\h^}[vW]* + [A - + ifcg^ ) 



+ [^ 1 -/ l W] M W(0), (4.5) 

where gw) is the function (| 1 . 14|) and we dropped, for the sake of brevity, the 
arguments (y,rj). Note that in writing down (|4. 5|) we used the fact that 



[(d v u)) 1 d t d v(j 



<P u 'i \(t,x;y,n)=(0,y,y,ri) 



compare with formula (|3.20p . 

Substituting formulae (|4.4p and (|4.5p into (|4.3|l we get 



where 



- h^)pWu { H(0) + P {l) B ( o 0) = 0, I ? j, 



(4.6) 



B, 



A ^ - lh% a + l (A 1 ) ya . na ) P« - ih<£P$ + i^O^pW. (4.7) 



The subscript in Pq^ indicates the degree of homogeneity in rj. 
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4.2 Part 2 of the proof of formula |L2| 

Our task in this subsection is to solve the linear algebraic system (|4.6p . (|4.2p for 
the unknowns u^\(0). 

It is easy to see that the unique solution to the system (|4.6[) , (|4.2p is 

Summation in (|4.8[) is carried out over all I different from j. 

4.3 Part 3 of the proof of formula ( 1TT201) 

Our task in this subsection is to calculate [£A-"(0)] su b- 
We have 

[^0)W=t£l(0)-~Pft a . (4.9) 

Here the sign in front of | is opposite to that in (| 1 . 1 5|) because the way we write 
U^{0) is using the dual symbol. 

Substituting ((4781) and (|4~71) into (|4~9| we get 

[c/ w (o)] sub =-^« Q +5: 1 



2 y "<> ^ /i0') - «0 



+ P«p + i(ii)v^)-P (i) - 42^- + W^]) 
v P(')A sub P^)+pO)A sub p(') w (J) v gg N 

^ - WO 2 V + 2^ h (j) - wo J ' ^- iUj 



where 



+ p 0) [(A W ^ (0 - 2<o p(0 + 2(^0^ P« ] . 
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We have 

Gji = 2P {1) {A 1 ,P { J ) } + 2P {J) {A 1 ,P {1) } 

+ P®[(Ai - h^) yar]a pW + 2(Ai - fcW)) fa pW] 
+ pWp! - h®) y « Va P® + 2(Aj - 
= 2PW { Ax , P« } + 2PV) { A x , pW } - P« { Ai - /i w , P w } - P {J) { A x - feW , P W } 
+ pOp! - ^)) ya , a P« + (A, - hW) Va P$ + (A, - ZiW^pW] 
+ P«[(Ai - ftW) w «, a PW + (Ax - ftW^pW + (Ax - ^),,aPW] 
= P®{A l + h®,pW} + pU){A 1 + h®,P w } 
- pW(Ax - ^ W )P^,„ - P®(Ai - fe«)Pft a 
= pW { Ax + fcW , P« } + P« { Ax + fcW , P« } 
- pW^CQ - fcW)P$,„ - /'(•') (//^ - /i^P^ 
= PW {Ax , P« }+P« { , P«}+(/i« - ) (P ( °P^ -P W P$ Va ) , 
so formula (|4. 10p can be rewritten as 

[t^ww = i(-p^„ +^(p«p« )7a -p«p,«j) 

, 1 y, PW (2A sub pO') + + ftW , }) + pO') (2^ sub P« + i{ A x + fc« , P« }) 



2 - fr« 



But 



(4.11) 



= (/ - pW)p^ o - p«(/ - p w w - p^, 

so formula (|4.1ip can be simplified to read 
[£/ (j) (0)]sub 

PW (2A sub P« + i { Ax + fcO') , p(i) }) + pU) (2A suh P^ +t{A 1 + h® , P w }) 

9 ' 



2 ^ /iCj) - fc(Q 



4.4 Part 4 of the proof of formula (11.201) 

Our task in this subsection is to calculate tr[{7^^(0)] su b- 



(4.12) 
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Formula (|4.12[) implies 

^«.^'E f " l4f 2^' 1 ^' 1 (-) 

Put Ai = 2fc hWpW and observe that 

• terms with the derivatives of h vanish and 

• the only k which may give nonzero contributions are k — j and k = I. 
Thus, formula (|4.13[) becomes 

x(h^ [pW {PW) , pO') } + p(i) {pCO ) p(0 }] [pCO {p(0 , p(i) }+p0') {p« , pW }] ) . 

(4.14) 

We claim that 

tr (p(0 {pCO ; pO) }) = tr (p(i) {pO) j p(0 }) 

= -tr(pW{pW,P«}) = -tr(P w {P (i) ,P (i) }) = 

These facts are established by writing the orthogonal projections in terms of the 
eigenvectors and using, if required, the identities 

+ [v w ]*v$ = 0, [v®]*v& + [v®]*v$ = 0, 

[v$]*v® + [ V U)]* V ® = 0, [v$]*vW + [v®]*v® = 0. 
In view of the identities (|4.15l) formula (14. 14[) can be rewritten as 

tr[[/«(0)] sub = itr^P«{P« P«} 

= i tr({P« , P«) } - pO) {pO') , P« }) = -i tr(P« {P^ , pW }) . (4.16) 
It remains only to simplify the expression in the RHS of (|4.16[) . We have 
tr (pW){p(j) } pCj)j.) = {[«(j)]* jt ,(j)j 

+ [([« tf) ]*«g?)([t;W]*t;W) - ([^l^JC^]*^)] 

= {[ V U)]*, V U)}. (4.17) 

Formulae (|4TT^|) and (|4~P7|) imply formula (fTT^D)) . 
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5 U(l) connection 



In the preceding Sections we presented technical details of the construction of 
the propagator. We saw that the eigenvectors of the principal symbol, v^(x, £), 
play a major role in this construction. As pointed out in Section [1] each of these 
eigenvectors is defined up to a U(l) gauge transformation (j 1 . 1 8(1 . 19|) . In the 
end, the full symbols (jl.lip of our oscillatory integrals U^'{t) do not depend on 
the choice of gauge for the eigenvectors ti* (x, £). However, the effect of the gauge 
transformation (|1.18p . (| 1 . 1 0|) is not as trivial as it may appear at first sight. We will 
demonstrate in this section that the gauge transformation (|1.18p . (|1.19p shows up, 
in the form of invariantly defined curvature, in the lower order terms u_\(t', y, rf) 
of the symbols of our oscillatory integrals U^\t). More precisely, we will show 
that the RHS of formula (|1.20l) is the scalar curvature of a connection associated 
with the gauge transformation (|1.18p . (|1.19p . Further on in this section, until the 
very last paragraph, the index j enumerating eigenvalues and eigenvectors of the 
principal symbol is assumed to be fixed. 

Consider a smooth curve T C T'M connecting points (y, rf) and (x, £). We write 
this curve in parametric form as (z(t),£(t)), t £ [0, 1], so that (z(0),£(0)) = {y,rf) 
and 0(1), C(l)) = (x,0- Put 

w(t) :=e^v^(z(t),((t)), (5.1) 

where </> : [0, 1] — > M is an unknown function which is to be determined from the 
condition 

iw*w = (5.2) 

with the dot indicating the derivative with respect to the parameter t. Substituting 
(|5.ip into (|5.2p we get an ordinary differential equation for which is easily solved, 
giving 

0(1) = 0(0) + f\z a (t) P a (z(t),{(t)) + C 7 (t) CT(z{t), C(t))) dt 
Jo 

= 0(0) + J (P a dz a + Q^dCj) , (5.3) 

where 

P a :=i[vW]*v$!, CP :=i[v®]*vg. (5.4) 

Note that the 2n-component real quantity (P a , Q 7 ) is a covector field (1-form) on 
T'M. This quantity already appeared in Section [T] as formula (|1.2ip . 
Put f(y,ri) := e*'^ \ f{x,£) := e 1 ^ and rewrite formula (gSD as 

f{x,0 = /Me^^+Q^f-r). (5.5) 

Let us identify the group U(l) with the unit circle in the complex plane, i.e. with 
/ G C, |/| = 1. We see that formulae (|5.5p and (|5.4p give us a rule for the parallel 
transport of elements of the group U(l) along curves in T'M . This is the natural 
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U(l) connection generated by the normalised field of columns of complex- valued 
scalars 

v^(z,0 = {v{ j \z,0 ... vtt(zX)) T - (5.6) 

Recall that the T appearing in formula (|5.5p is a curve connecting points (y, rf) 
and (#,£), whereas the v^'(z, £) appearing in formulae (|5.4p and (|5 . 6[) enters our 
construction as an eigenvector of the principal symbol of our mxm matrix pseudo- 
differential operator A. 

In practice, dealing with a connection is not as convenient as dealing with the 
covariant derivative V. The covariant derivative corresponding to the connection 
(|5.5|l is determined as follows. Let us view the (x, £) appearing in formula ()5.5[) as 
a variable which takes values close to (y, rf), and suppose that the curve T is a short 
straight (in local coordinates) line segment connecting the point (y, rf) with the 
point (x, £). We want the covariant derivative of our function f(x, £), evaluated at 
(y,Tj), to be zero. Examination of formula (15 . 5[) shows that the unique covariant 
derivative satisfying this condition is 

V a :=d/dx a -iP a (x,0, V* -.= 8/8$, -iCP(x,t). (5.7) 

We define the curvature of our U(l) connection as 

VaV/3-V/jVa V a V 5 -X7 s X7 c 



n ._ _ ■ v a » p vpv Q v Q v v v qj 1 /cr o\ 

* V 7 V,3 - V3V 7 V 7 V 5 -V 5 V 7 '' 



It may seem that the entries of the (2n) x (2n) matrix (|5.8|) are differential opera- 
tors. They are, in fact, operators of multiplication by "scalar functions". Namely, 
the more explicit form of (|5.8|) is 





_ dPp_ 


dP a _ dQ s 




dx a 


d£s dx a 


dQ~" 


dPff 


dQ 1 _ dQ s 


dx& 




d£s d£ 7 



E> _ dxP dx a d£s dx" \ /r q\ 

11 — do~< 9P« do-i do 6 ■ 



The (2n) x (2n) - component real quantity (|5.9p is a rank 2 covariant antisymmetric 
tensor (2-form) on T'M. It is an analogue of the electromagnetic tensor. 

Substituting (|5.4I) into (|5.9p we get an expression for curvature in terms of the 
eigenvector of the principal symbol 

Examination of formula (|5.10|) shows that, as expected, curvature is invariant 
under the gauge transformation (j!.18[) . (|1.19D . 

It is natural to take the trace of the upper right block in (|5.8p which, in the 
notation (|1.16p . gives us 

-i(V Q V Q - V a V Q ) = -i{[v^]*,v (j) }. (5.11) 

Thus, we have shown that the RHS of formula f| 1 . 20[) is the scalar curvature of our 
U(l) connection. 
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We end this section by proving, as promised in Section[TJ formula (| 1 . 22[) without 
referring to microlocal analysis. In the following arguments we use our standard 
notation for the orthogonal projections onto the eigenspaces of the principal sym- 
bol, i.e. we write pW := v {k ^[v^}*. We have tr{P<>'),pW} = and £, P (i) = 1 
which implies 

= ^tr(p( () {P«,p( J }}) 

hi 

= J2 t *( pU) {P U \P U} })+ J2 te(P {l) {P U \P U} })- (5-12) 

3 hj- l^j 

But, according to formula f|4 . 1 5[) . for I ^ j we have 

tr (p(0{pW ;P W}) = -tr(pW{pW,P«}), 

so the expression in the last sum in the RHS of (I5.12j) is antisymmetric in the 
indices l,j, which implies that this sum is zero. Hence, formula (|5.12[) can be 
rewritten as J^tr (Pw{pW,pW}) = 0. It remains only to note that, according 

3 

to formula ijCTF]) . tr(P^{pW, P^}) = {[v^]*,v^}. 

6 Singularity of the propagator at t = 

Following the notation of [TB], we denote by 

•?W/(A)] = f(t) = J e- ux f(X)d\ 
the one-dimensional Fourier transform and by 

Ff-lxlM = fW = W 1 J e ltX f(t)dt 

its inverse. 

Suppose that we have a Hamiltonian trajectory (xtO (t; y, 77), (t; y, 77)) and 
a real number T > such that x^ (T; y, 77) — y. We will say in this case that we 
have a loop of length T originating from the point y £ M. 

Remark 6.1 There is no need to consider loops of negative length T because, 
given a T > 0, we have x^(T;y,ri + ) = y for some r/ + £ TyM if an d only if we 
have x^ (— T; y, 77^) = y for some r)~ £ TyM. Indeed, it suffices to relate the rf 1 
in accordance with if* 1 = ^-^(iT;?/,^). 

Denote by T"' C K the set of lengths T > of all possible loops generated 
by the Hamiltonian ftW. Here "all possible" refers to all possible starting points 
(y,n) £ T'M of Hamiltonian trajectories. It is easy to see that ^ ■ We put 




infT w if T (i) ^0, 
+00 if T (i) = 0. 
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In the Riemannian case (i.e. the case when the Hamiltonian is a square root of 
a quadratic polynomial in £) it is known [HI H"2"] that there is a loop originating 
from every point of the manifold M and, moreover, there is an explicit estimate 
from above for the number . We are not aware of similar results for general 
Hamiltonians. 

We also define T := min T^'). 

j — l,...,m+ 

Remark 6.2 Note that negative eigenvalues of the principal symbol, i.e. Hamil- 
tonians hU'(x,£) with negative index j = — 1, . . . , —mT, do not affect the asymp- 
totic formulae we are about to derive. This is because we are dealing with the case 
A — > +00 rather than A — > — 00. 

Denote by 

u(t,x,y) :=J2 e ~ itXkv k( x )l v k(y)}* (6-1) 
k 

the integral kernel of the propagator . The quantity (|6.ip can be understood 
as a distribution in the variable t G R depending on the parameters x,y G M. 
The main result of this section is the following 

Lemma 6.1 Let p : K — > C be an infinitely smooth function such that 

suppp C (— T, T), (6.2) 

p(0) = 1, (6.3) 
p'(0) = 0. (6.4) 
Then, uniformly over y G M , we have 

F~_l x [p(t) tr u(t, y, y)] = n a(y) A" -1 + (n - 1) 6(y) A"~ 2 + 0(A"~ 3 ) (6.5) 



as A — ^ +00. 27ie densities a(y) and b(y) appearing in the RHS of formula i6.5]) 
are defined in accordance with formulae $1.23\) and \1.2J$ . 

Proof Denote by (S*M)^ the (n— l)-dimensional unit cosphere in the cotan- 
gent fibre defined by the equation h^'(y,rf) = 1 and denote by d(S*M)^' the 
surface area element on (S*M)^' defined by the condition dr\ = d{S*M)^ dh^\ 
The latter means that we use spherical coordinates in the cotangent fibre with the 
Hamiltonian ftW) playing the role of the radial coordinate, see subsection 1.1.10 of 
[15] for details. In particular, as explained in subsection 1.1.10 of [18], our surface 
area element d{S*M)^ is expressed via the Euclidean surface area element as 

/ ™ \-i/2 
d(S*M)^ = ^(^(l/,*?)) x Euclidean surface area element . 

Denote also d(S*M)^ := (2%)- n d(5*M)W . 
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According to Corollary 4.1.5 from [T5] we have uniformly over y £ M 
•?vL[p(*) tru (*>2/>y)] = 

m+ 

(c 0) (y) A"" 1 + (y) A"~ 2 + e« (y) A"~ 2 ) + 0(A"- 3 ) , (6.6) 



where 



C W(y)= j tr (0; y, V ) d(S* v M)^ , (6.7) 

(ri- 1) | tr^z4 j) (0;y,r ? ) + ^{4 j) | t=0 ,/ l «}(y,r ? )^(5;M)W, (6.8) 

(S;Af)0) 

e (j) (2/)= / tr[trW(0)] Bub (i/,»?)d(5*A0 W) . (6.9) 

(S*Af)U) 

Here uffi (t; y, rf) is the principal symbol of the oscillatory integral (|2 . 1 1) and ii^ (t; y, rf) 
is its time derivative. Note that in writing the term with the Poisson bracket in 
(|6.8p we took account of the fact that Poisson brackets in |18| and in the current 
paper have opposite signs. 

Observe that the integrands in formulae (|6.7[) and (|6.8[) are positively homo- 
geneous in -q of degree 0, whereas the integrand in formula (|6.9[) is positively 
homogeneous in r\ of degree —1. In order to have the same degree of homogeneity, 
we rewrite formula (|6.9p in equivalent form 



e (j) (2/)= / (h (j hY[U^(0)} suh )(y,f 1 )d(S;M)^ . (6.10) 

Switching from surface integrals to volume integrals with the help of formula 
(1.1.15) from [TB], we rewrite formulae ([6T7| . (|^5j) and (pTTUj) as 

c u) (y)=n J trutf ) (Q;y,r))dv, (6.11) 

d (j \y) = n(n-l)x 

j tr(-iutt\a;y,r,) + ~{uiX =0 ,h^(y,v)^<IV, (6-12) 

/iU)(y,77)<l 
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e^(y) = n I (h® tr[L/«(0)] sub ) (y, rj) dr) . (6.13) 

h(J) (y,ri)<l 

Substituting formulae (j!.12|l and (|1.14[) into formulae (|6.1ip and (|6.12ll we get 




(6.14) 



d^(y) = - n (n-l)x 

[v^]*A suh v^ - -{[v^F^-h®^®}) {y,n)dn. (6.15) 

Substituting formula (|1.20|) into formula (|6.13| we get 

e {l) (y) = -ni J (h^{[v u) }*, v u} })(y, rj) tir) . (6.16) 

W>(y,r))<l 

Substituting formulae (|6.14[) - (|6.16p into formula (|6.6p we arrive at (|6.5[) . □ 

Remark 6.3 The proof of Lemma \6.1\ given above was based on the use of Corol- 
lary 4-1-5 from JTBjj. In the actual statement of Corollary 4. 1.5 in f!8f uniformity 
in y £ M was not mentioned because the authors were dealing with a manifold 
with a boundary. Uniformity reappeared in the subsequent Theorem 4-2.1 which 
involved pseudodifferential cut-offs separating the point y from the boundary. 

7 Mollified spectral asymptotics 

Theorem 7.1 Let p : M — > C be a function from Schwartz space S(M) whose 
Fourier transform p satisfies conditions Ii6.^) - [6^4\ l. Then, uniformly over x £E M , 
we have 

fe(\-v,x,x)p(p)dv = a(x)\ n +b(x)\ n - 1 + r^ n ~ 2) * " " J' (7.1) 
J ^O(lnA) if n = 2, 

as A — > +oo. The densities a(x) and b(x) appearing in the RHS of formula i7.1]) 
are defined in accordance with formulae il.23\) and \1.2J$ . 

Proof Our spectral function e(A, x, x) was initially defined only for A > 0, see 
formula (|1.4j) . We extend the definition to the whole real line by setting 

e(A, x, x) := for A < 0. 
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Denote by e'(A, x, x) the derivative, with respect to the spectral parameter, of 
the spectral function. Here "derivative" is understood in the sense of distributions. 
The explicit formula for e'(A, x, x) is 



+00 



e'(A,x,x) :=^2\\v k (x)fS(X-X k ), (7.2) 



k=l 

Formula (|7.2[) gives us 

+00 



e'(A - p,x,x) p{p)dp = ^ \\vk{x)\\ 2 p{\ - Afe). (7.3) 
fe=i 

Formula (|7.3p implies, in particular, that, uniformly over x £ M , we have 

e'(A — /1, £E, x) dp — 0(|Ap°°) as A — >• — 00 , (7-4) 



where 0(|A| 00 ) is shorthand for "tends to zero faster than any given inverse power 
of |A|". 

Formula (|7.3|) can also be rewritten as 

e'(X-fi,x,x)p(fi)dfj, = J^ x [p(t)tru(t,x,x)] ||u fc (x)|| 2 p(A - A fc ) , (7.5) 

fc<0 

where the distribution u(t,x,y) is defined in accordance with formula (|6.1I) . Clearly, 
we have 

^|| Ufe (x)[| 2 /9 (A-A fe )=0(A- 00 ) as A -> +00 . (7.6) 

fe<0 

Formulae (|7.5p . (I7.6[) and Lemma IUTT1 imply that, uniformly over x 6 M, we have 

e (A — /i, a;, a;) p(p) dp, = 

na(a;)A™- 1 + (n-l)fo(x)A™- 2 +0(A™- 3 ) as A -> +00 . (7.7) 
It remains to note that 

— [ e(X- p,x,x) p(p)dp = ( e'{X- p,x,x) p(p)dp. (7.8) 



dA . 

Formulae dTTHJ), {HJ) and ([721) imply flU]). □ 

Theorem 7.2 Lei p : 1 -> C 6e a function from Schwartz space 5(R) whose 
Fourier transform p satisfies conditions i6.2]) - [K^\ ). Then we have 

N(X-p)p(p)dp = aX n + bX n - 1 + \° {Xn ^ ) lf n " 3 ' (7.9) 



O(lnA) if n = 2, 

as A — > +00. The constants a and b appearing in the RHS of formula \7.9\) are 
defined in accordance with formulae hi. 8)) , hi. 23]) , lil.9\) and 11.24]) - 
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Proof Formula (|7.9p follows from formula (|7.ip by integration over M, see also 
formula ([L5]) . □ 

In stating Theorems 17.11 and 17.21 we assumed the molliher p to be complex- 
valued. This was done for the sake of generality but may seem unnatural when 
mollifying real- valued functions e(A, x, x) and N(X). One can make our construc- 
tion look more natural by dealing only with real-valued mollifiers p. Note that 
if the function p is real-valued and even then its Fourier transform p is also real- 
valued and even and, moreover, condition (|6.4[) is automatically satisfied. 

8 Unmollified spectral asymptotics 

In this section we derive asymptotic formulae for the spectral function e(X,x,x) 
and the counting function N(X) without mollification. The section is split into 
two subsections: in the first we derive one-term asymptotic formulae and in the 
second — two-term asymptotic formulae. 

8.1 One-term spectral asymptotics 
Theorem 8.1 We have, uniformly over x G M , 

e(\,x,x) = a(x)\ n + 0(X n - 1 ) (8.1) 

as X — > +oo. 

Proof The result in question is an immediate consequence of formulae (17. 8p , 
(|7.7p and Theorem 1 7 . 1 1 from the current paper and Corollary B.2.2 from [18]. □ 

Theorem 8.2 We have 

N(X) = a\ n + 0{X n ' 1 ) (8.2) 

as X ^ +oo. 

Proof Formula (|8.2p follows from formula (|8.ip by integration over M, see also 
formula HH). □ 

8.2 Two-term spectral asymptotics 

Up till now, in Section [7] and subsection 18. 1[ our logic was to derive asymptotic 
formulae for the spectral function e(A, x, x) first and then obtain corresponding 
asymptotic formulae for the counting function -/V(A) by integration over M . Such 
an approach will not work for two-term asymptotics because the geometric condi- 
tions required for the existence of two-term asymptotics of e(A, a;, x) and N(X) will 
be different: for e{X,x,x) the appropriate geometric conditions will be formulated 
in terms of loops, whereas for N(X) the appropriate geometric conditions will be 
formulated in terms of periodic trajectories. 

Hence, in this subsection we deal with the spectral function e(X,x,x) and the 
counting function N(X) separately. 
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In what follows the point y 6 M is assumed to be fixed. 

Denote by ILy the set of normalised (h^\y, rj) = 1) covectors 77 which serve as 
starting points for loops generated by the Hamiltonian hy>. Here "starting point" 
refers to the starting point of a Hamiltonian trajectory (x^'(t;y, T}) , £^ (t; y , rj)) 
moving forward in time (t > 0), see also Remark 16. II 

The reason we are not interested in large negative t is that the refined Fourier 
Tauberian theorem we will be applying, Theorem B.5.1 from \TE\, does not require 
information regarding large negative t. And the underlying reason for the latter 
is the fact that the function we are studying, e(A, x, x) (and, later, N(X)), is 
real-valued. The real-valuedness of the function e(\,x,x) implies that its Fourier 
transform, e(t, x, x), possesses the symmetry e(— t,x, x) = e{t,x,x). 

The set II^ is a subset of the (n — l)-dimensional unit cosphere (S*M)^ and 
the latter is equipped with a natural Lebesgue measure, see proof of Lemma |6. 11 
It is known, see Lemma 1.8.2 in [18j . that the set Tly ^ is measurable. 

Definition 8.1 A point y £ M is said to be nonfocal if for each j = 1, . . . ,m + 
the set Hy^ has measure zero. 

With regards to the range of the index j in Definition 18.11 as well as in sub- 
sequent Definitions I8.2H8.41 see Remark 16.21 

We call a loop of length T > absolutely focused if the function 

\x^(T;y,r,)-y\ 2 

has an infinite order zero in the variable 77, and we denote by (Hy)^ the set of 
normalised (h^(y, rf) = 1) covectors 77 which serve as starting points for absolutely 
focused loops generated by the Hamiltonian hV>. It is known, see Lemma 1.8.3 in 
|18j . that the set (Hy)^ is measurable and, moreover, the set U y 3) \ (II£)W has 
measure zero. This allows us to reformulate Definition 18.11 as follows. 

Definition 8.2 A point y 6 M is said to be nonfocal if for each j = 1, ... ,rn + 
the set (II")' 3 ' has measure zero. 

In practical applications it is easier to work with Definition 18.21 because the set 
(n°)( J ) is usually much thinner than the set liy . 

In order to derive a two-term asymptotic formula for the spectral function 
e(A, x, x) we need the following lemma (compare with Lemma |6.1[) . 

Lemma 8.1 Suppose that the point y 6 M is nonfocal. Then for any complex- 
valued function 7 G C^°(R) with supply C (0, +00) we have 

T-^ x mt)tru(t,y,y)} = o(A"- 1 ) (8.3) 

as A — > +00. 

Proof The result in question is a special case of Theorem 4.4.9 from [T5]. □ 
The following theorem is our main result regarding the spectral function e(A, x, x). 
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Theorem 8.3 If the point x £ M is nonfocal then the spectral function e(A, x, x) 
admits the two-term asymptotic expansion (1.6\) as X — > +00. 

Proof The result in question is an immediate consequence of formulae (17. 7p , 
Theorem 17.11 and Lemma 18.11 from the current paper and Theorem B.5.1 from 

HSI.n 

We now deal with the counting function N(X). 

Suppose that we have a Hamiltonian trajectory (x^(t;y ) rf),^'(t' } y ) rf)) and 
a real number T > such that (x^> (T; y, rj),^' (T; y, rf)) — (y, vf). We will say in 
this case that we have a T-periodic trajectory originating from the point (y, 77) G 
T'M. 

Denote by {S*M)V> the unit cosphere bundle, i.e. the (2n — l)-dimensional 
surface in the cotangent bundle defined by the equation h^(y,if) = 1. The 
unit cosphere bundle is equipped with a natural Lebesgue measure: the (2n — 1)- 
dimensional surface area element on (S*M)^ is dy d(S*M)^ where d{S*M)^ 
is the (n — l)-dimensional surface area element on the unit cosphere (S*M)^', see 
proof of Lemma 16.11 

Denote by 11^ the set of points in (S*M)^ which serve as starting points for 
periodic trajectories generated by the Hamiltonian h^K It is known, see Lemma 
1.3.4 in [18], that the set II^) is measurable. 

Definition 8.3 We say that the nonperiodicity condition is fulfilled if for each 
j = 1, . . . , m + the set has measure zero. 

We call a T-periodic trajectory absolutely periodic if the function 

\xU)(T;y,T})-y\ 3 + \zW(T ]y ,Ti)-ri\ 3 

has an infinite order zero in the variables (j/,77), and we denote by (II a )W the 
set of points in (S*M)V' which serve as starting points for absolutely periodic 
trajectories generated by the Hamiltonian h^K It is known, see Corollary 1.3.6 in 
|18j . that the set (H a )^ is measurable and, moreover, the set \ (n a )(j) has 
measure zero. This allows us to reformulate Definition 18.31 as follows. 

Definition 8.4 We say that the nonperiodicity condition is fulfilled if for each 
j = 1, . . . ,m + the set (JJ a )^ has measure zero. 

In practical applications it is easier to work with Definition 18.41 because the set 
(n°)w is usually much thinner than the set n^). 

In order to derive a two-term asymptotic formula for the counting function 
N(X) we need the following lemma. 

Lemma 8.2 Suppose that the nonperiodicity condition is fulfilled. Then for any 
complex-valued function 7 6 C^°(K) with supply C (0,+oc) we have 

I T^ x m)tTu(t,y 1 y)]dy = o(X n - 1 ) (8.4) 
Jm 

as X — > +00. 
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Proof The result in question is a special case of Theorem 4.4.1 from [T5]. □ 
The following theorem is our main result regarding the counting function N(X). 

Theorem 8.4 If the nonperiodicity condition is fulfilled then the counting func- 
tion N(X) admits the two-term asymptotic expansion il.7\) as A — >• +00. 

Proof The result in question is an immediate consequence of formulae (jl.5p , 
(|7.7p . Theorem 17.11 and Lemma [8.21 from the current paper and Theorem B.5.1 
from [18J. □ 



9 U(m) invariance 

We prove in this section that the RHS of formula (|1.24[) is invariant under unitary 
transformations (|1.25p , (|1.26p of our operator A. The arguments presented in this 
section bear some similarity to those from Section^ the main difference being that 
the unitary matrix-function in question is now a function on the base manifold M 
rather than on T'M. 

Fix a point x G M and an index j (index enumerating the eigenvalues and 
eigenvectors of the principal symbol) and consider the expression 

[v^}*A suh v^ 
fc«>(x,£)<l 

- U[v&]*,Ai - h {j) , } + —t—h^ { [vW]*,v& } W d£ , (9.1) 

2 ' n — 1 / 

compare with (|1.24p . We will show that this expression is invariant under the 
transformation ([OS)) , ([06) , 

The transformation (|1.25p , (|1.26p induces the following transformation of the 
principal and subprincipal symbols of the operator A: 

A t 1-4 RA X R*, (9.2) 

A suh ^ RA suh R* + l - {R xa (A 1 )^ a R* - R(Ai)^ a R* a ) . (9.3) 

The eigenvalues of the principal symbol remain unchanged, whereas the eigen- 
vectors transform as 

v U)^ Rv U)_ (9.4) 
Substituting formulae (|9.2l) - (|9.4p into the RHS of (|9.ip we conclude that the in- 
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crement of the expression (|9.1I) is 

J ^L [v U) r{R * Rxa{Al)ia { Al ) ia Rl a R)v^ 

- \ [[v {0) ]*R* xa R{A 1 - h.W)vg - - h^)R*R x a,v^ 

which can be rewritten as 

([v U) ]*R*<* RvU) ~ [v u) }*R*R x «v^ 

fiO*)( x ,f)<l 

- -^—^h { ^ {[v^]*R* xa Rv^ - [v^]*R*R x o,v^yj(x,0 d£ . 

In view of the identity R*R = I the above expression can be further simplified, so 
that it reads now 

h^[v (J) ]*R*R x o,v (J) 

- ^jh^ ([vW]*R*R x c, v V] + [v^]*R*R xa v^U(x,() d£ . (9.5) 

Denote B a (x) :— —iR*R x a and observe that this set of matrices, enumerated 
by the tensor index a running through the values 1, . . . ,n, is Hermitian. Denote 
also b a (x, £ ) :— [v^']*B a v^' and observe that these b a are positively homogeneous 
in £ of degree 0. Then the expression (j9.5[) can be rewritten as 

- / (^-^ (I, 1:)<^- 

Lemma 4.1.4 and formula (1.1.15) from [18] tell us that this expression is zero. 

10 Spectral asymmetry 

In this section we deal with the special case when the operator A is differential 
(as opposed to pseudodifferential). Our aim is to examine what happens when we 
change the sign of the operator. In other words, we compare the original operator 
A with the operator A := —A. In theoretical physics the transformation A n> — A 
would be interpreted as time reversal, see equation (|1.3[) . 

It is easy to see that for a differential operator the number m (number of 
equations in our system) has to be even and that the principal symbol has to have 





35 



the same number of positive and negative eigenvalues. In the notation of Section [T] 
this fact can be expressed as to = 2m + = 2mT . 

It is also easy to see that the principal symbols of the two operators, A and A, 
and the eigenvalues and eigenvectors of the principal symbols are related as 



Ai(x,0 


= A 1 {x,-0, 


(10.1) 






(10.2) 




= v( j) {x,-£), 


(10.3) 



whereas the subprincipal symbols are related as 

A mh (x) = -A suh (x). (10.4) 
Formulae (fl~2"3l) , ([T23|) . (TLlTj) , (fl~T5)) and (TilTTjl - pIHl) imply 

a(x) = a{x), b(x) = -b(x). (10.5) 
Substituting (IT03)) into (TT8]) and (TO]) we get 

a = a, b=-b. (10.6) 

Formulae (|1.7p and (|10.6p imply that the spectrum of a generic first order 
differential operator is asymmetric about A = 0. This phenomenon is known as 
spectral asymmetry [TJ [5J [3J 0] . 

If we square our operator A and consider the spectral problem A 2 v = X 2 v, then 
the terms ±&A n_1 cancel out and the second asymptotic coefficient of the counting 
function (as well as the spectral function) of the operator A 2 turns to zero. This 
is in agreement with the known fact that for an even order semi-bounded matrix 
differential operator acting on a manifold without boundary the second asymptotic 
coefficient of the counting function is zero, see Section 6 of [20] and |15j . 

11 Bibliographic review 

To our knowledge, the first publication on the subject of two-term spectral asymp- 
totics for systems was Ivrii's 1980 paper [5] in Section 2 of which the author stated, 
without proof, a formula for the second asymptotic coefficient of the counting func- 
tion. In a subsequent 1982 paper [7] Ivrii acknowledged that the formula from [6] 
was incorrect and gave a new formula, labelled (0.6), followed by a "proof". In 
his 1984 Springer Lecture Notes [8] Ivrii acknowledged on page 226 that both his 
previous formulae for the second asymptotic coefficient were incorrect and stated, 
without proof, yet another formula. 

Roughly at the same time Rozenblyum [13] also stated a formula for the second 
asymptotic coefficient of the counting function of a first order system. 

The formulae from [5], [7] and [T3J are fundamentally flawed because they are 
proportional to the subprincipal symbol. As our formulae (|1.9p and (|1.24p show, 
the second asymptotic coefficient of the counting function may be nonzero even 
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when the subprincipal symbol is zero. This illustrates, yet again, the difference 
between scalar operators and systems. 

The formula on page 226 of [8| gives an algorithm for the calculation of the 
correction term designed to take account of the effect described in the previous 
paragraph. This algorithm requires the evaluation of a limit of a complicated 
expression involving the integral, over the cotangent bundle, of the trace of the 
symbol of the resolvent of the operator A constructed by means of pseudodiffcr- 
ential calculus. This algorithm was revisited in Ivrii's 1998 book, see formulae 
(4.3.39) and (4.2.25) in 0. 

The next contributor to the subject was Safarov who, in his 1989 DSc The- 
sis JT5], wrote down a formula for the second asymptotic coefficient of the count- 
ing function which was "almost" correct. This formula appears in [TB] as for- 
mula (2.4). As explained in Section [1] Safarov lost only the curvature terms 
I h^{[v Vi]*,v V'}-. Safarov's DSc Thesis [T^] provides arguments which are 
sufficiently detailed and we were able to identify the precise point (page 163) at 
which the mistake occurred. 

In 1998 Nicoll rederived [llj Safarov's formula (|1.12p for the principal symbols 
of the propagator, using a method slightly different from [16], but stopped short 
of calculating the second asymptotic coefficient of the counting function. 

In 2007 Kamotski and Ruzhansky [10] performed an analysis of the propa- 
gator of a first order elliptic system based on the approach of Rozenblyum [15] . 
but stopped short of calculating the second asymptotic coefficient of the counting 
function. 

One of the authors of this paper, Vassiliev, considered systems in Section 6 
of his 1984 paper [20]. However, that paper dealt with systems of a very special 
type: differential (as opposed to pseudodifferential) and of even (as opposed to 
odd) order. In this case the second asymptotic coefficients of the counting function 
and the spectral function vanish, provided the manifold does not have a boundary. 
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